Remark on the irrationality of the Brun's constant 



Marek Wolf 



e-mail:mwolf@ift. uni.wroc.pl 



Abstract 



We have calculated numerically geometrical means of the denominators 
of the continued fraction approximations to the Brun constant 82- We 
get values close to the Khinchin's constant. Next we calculated the n-th 
square roots of the denominators of the n-th convergents of these con- 
tinued fractions obtaining values close to the Khinchin-L'evy constant. 
These two results suggests that £>2 is irrational, supporting the common 
believe that there is an infinity of twins. 



Very well known open problem in number theory is the question whether there 
exist infinitely many twin primes p, p + 2. In 1919 Brun [3] has shown that the sum 
of the reciprocals of all twin primes is finite: 




thus leaving the problem not decided. Sometimes 5 is included in (II]) only once, but 
here we will adopt the above convention. The sum ([lj is called the Brun constant 



Let ^(x) denote the number of twin primes smaller than x. Then the conjecture 
B of Hardy and Littlewood [5] on the number of prime pairs p,p + d applied to the 
case d = 2 gives, that 



[lOj. 




where C 2 is called "twin constant" and is defined by the following infinite product: 
C 2 = 2T](l-- — = 1.32032363169... (3) 

There is a large evidence both analytical and experimental in favor of ^2§. Besides 
the original circle method used by Hardy and Littlewood [5] there appeared the 
paper [8] where another heuristic arguments were presented. In May 2004, in a 
preprint publication [lj Arenstorf attempted to prove that there are infinitely many 
twins. Arenstorf tried to continue analytically to 3?(s) = 1 the difference: 

T(s) - C 2 /(s - 1) (4) 

where the function 

T(s) = J2 A ( n ~ l)A(n+ l)n~ s (3ft(s) > 1). (5) 

n>3 

However shortly after an error in the proof was pointed out by Tenenbaum [TT] . For 
recent progress in the direction of the proof of the infinite number of twins see [6] . 
Because there is no doubt that twins prime conjecture is true the Brun's constant 
should be irrational. 

The series (JTJ is very slowly convergent and there is a method based on the ^ 
to extrapolate finite size approximations 

«■<*>= £ H + < 6 > 

q,q + 2 both prime - 1 ± 
q<x 

to infinity [2]: 

Bt\x)=B 2 (x) + —^- (7) 
log(x) 

In this way from the straight sieving of primes up to x = 3 x 10 15 T. Nicely [7] gives 
fi 2 (oo) (3 x 10 15 ) = 1.9021605823 ± 0.0000000008 (8) 
while P. Sebah [9] from computer search up to 10 16 gives 

# 2 M(10 16 ) = 1.902160583104. (9) 

If there is an infinity of twins, as the formula ^ asserts, then the Brun's constant 
should be an irrational number. Vice versa if the Brun's constant is irrational then 
there is an infinity of twins. 

There exists a method based on the continued fraction expansion which allows 
to detect whether a given number r can be irrational or not. Let 

r = [a Q ;a 1 ,a 2 ,a 3 , . . .] = a H (10) 



a\ H 

a>2 



a 3 + 



2 



be the continued fraction expansion of the real number r, where ao is an integer and 
all cik, A; = 1,2,... are positive integers. Khinchin has proved that 

j 00 f 1 ^ log 2 m 

Urn (a 1 a 2 ...a n )" = T] + — -} = K « 2.685452001 .. . (11) 

n^oo ' J-x m(m + 2) J 

is a constant for almost all real r, see e.g. [H §1-8]. The exceptions are rational 
numbers, quadratic irrationals and some irrational numbers too, like for example the 
Euler constant e = 2.7182818285. . ., but this set of exceptions is of the Lebesgue 
measure zero. The constant Kq is called the Khinchin constant. If the quantities 

K(n) = {a ia2 ...a n y (12) 

for r given with accuracy of some number of digits are close to Kq we can regard it 
as a hint that r is irrational. 

For the numerical value of B 2 given by ^ we get continued fraction containing 
23 terms: 



B 2 w 1.902160583104 = [1; 1, 9, 4, 1, 1, 8, 3, 4, 4, 2, 2, 2, 1, 35, 1, 1, 1, 2, 4, 4, 1, 2] 

1 + — (13) 



1 + 



1 



9 + 



1 + 



+ 



4+ 

We have calculated the geometrical means K(n) for the consecutive truncations of 
the continued fraction (13) for n = 7, 8,... 23. The results are presented in the 
second column of Table 1 and in Fig.l for n = 3, 4, . . . , 23 and K{n) are fluctuating 
around K , suggesting B 2 is indeed irrational. 
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TABLE I 



n 


K(n) 


Lin) 


7 


2.5697965868506505913 


3.038522491494629559394 


8 


2.6272534028385753462 


3.049597387327883630584 


9 


2.7689921001973389000 


3.185356988368901235392 


10 


2.8844991406148167646 


3.287467122665676043642 


11 


2.7807783406318186810 


3.163215603753400545370 


12 


2.6986960551881232673 


3.090263509588119664954 


13 


2.6321480259049849216 


3.026751930527273995360 


14 


2.4433109823349210157 


2.854823973874123136328 


15 


2.9549936641069114440 


3.419389588999534265980 


16 


2.7490720167650636040 


3.156107031770304366266 


17 


2.5806968087738011680 


3.064762600086287594909 


18 


2.4407155166256854482 


2.939414938043217059705 


19 


2.4138613420275263905 


2.923369271618110487047 


20 


2.4788881641188601981 


2.986079858472829061265 


21 


2.5389087172991914976 


3.038476735121596843751 


22 


2.4287240750956647922 


2.911134637211735156321 


23 


2.4073773469514141247 


2.848836285543623055400 



Let the rational p n /q n be the n-th partial convergent of the continued fraction: 

— = [a Q ;a 1 ,a 2 ,a 3 ,. . . ,a n ]. (14) 

For almost all real numbers r the denominators of the finite continued fraction 
approximations fulfill: 

lim {q n {r)) 1/n = e" 2/121n2 = L = 3.275822918721811 . . . (15) 

where L is called the Khinchin — Levy's constant [J, §1.8]. Again the set of excep- 
tions to the above limit is of the Lebesgue measure zero and it includes rational 
numbers, quadratic irrational etc. 

From (13) we get the following sequence of convergents p n (^2) I '<?n(£>2): 

1 2 19 78 97 175 1497 4666 20161 85310 190781 466872 
1' 1' 10' 41' 51' "92"' ~787~' 2453' 10599' 44849' 100297' 245443 ' 

1124525 1591397 56823420 58414817 115238237 173653054 



591183 ' 836626 ' 29873093' 30709719' 60582812 ' 91292531 ' 
462544345 2023830434 8557866081 10581696515 19139562596 \ 
243167874' 1063964027' 4499023982' 5562988009 ' 10062011991 J 
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From these denominators q n (B 2 ) we can calculate the quantities L(n): 

L{n) = (q n (B 2 )) 1/n , n = 1,2,..., 23 (16) 

The obtained values of L(n) for n > 7 are presented in the third column of the Table 
1 and are shown in the Fig.2. These values scatter around the red line representing 
the Khinchin — Levy's constant again suggesting that B 2 is irrational. 

In conclusion we can say that to draw firmer statement much more digits of the 
Brun's constant are needed. 



Acknowledgement: The calculation were performed using the free package 
PARI/GP [12] developed especially for number theoretical purposes and which has 
built in a set of function to deal with continued fractions. 
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Fig.l The plot of the consecutive geometrical means K(n), n = 3,4, . . . , 23. 
Although the number of available for the value ^ of ^points (N, K{n)) is rather 
moderate there are 5 sign changes of the difference K(n) — K . 
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Fig.2 The plot of the consecutive values of q n , n — 3, 4, . . . , 23. Although the 
number of available for the value ^ of ^points (n, L(n)) is rather moderate there 
are 6 sign changes of the difference L(n) — Lq. 
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